Theoretical study of NMR relaxation due to rattling phonons 
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Abstract 

We calculate the NMR relaxation rate due to quadrupolar coupling of the nucleus to a local, strongly anharmonic 
phonon mode. As a model potential for a "rattling" motion we consider a square-well potential. We calculate the free 
phonon Green's function analytically and derive the low and high temperature limits of the NMR relaxation rate. 
It is shown that the temperature dependence of the NMR relaxation rate possesses a peak in contrast to harmonic 
phonons but in qualitative agreement with a recent NMR study on KOS2O6. We discuss the influence of phonon 
renormalization due to electron-phonon interaction. 
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The new pyrochlore oxide superconductor 
KOS2O6 shows unusual behavior both in the nor- 
mal as well as in the superconducting state. Den- 
sity functional calculations have shown that the 
vibrations of the potassium ion are highly anhar- 
monic, allowing large excursions from its equilib- 
rium position within the surrounding Os-0 cage 
[1], consistent with recent X-ray studies [2]. Such 
a situation, in which a small ion can move anhar- 
monically within an oversized cage of surrounding 
atoms has been called "rattling" and also exists in 
other compounds [3,4,5]. The anharmonic dynamics 
of the potassium ion in KOS2O6 is currently being 
discussed as a possible origin for the anomalies seen 
in various experimental quantities, like e.g. specific 
heat [6], resistivity [7], thermal conductivity [8], 
and NMR relaxation [9]. 

Here we focus on recent NMR studies of KOS2O6. 
In particular, the nuclear spin-lattice relaxation rate 
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l/TiT at the potassium site can directly probe the 
unusual potassium dynamics through quadrupolar 
coupling of the K nucleus to the electric field gradi- 
ent, as has been shown recently [9] . The temperature 
dependence of the NMR relaxation rate 1/TiT was 
found to be highly anomalous, increasing at low tem- 
peratures, reaching a peak value around ~ 14 K and 
decreasing again at higher temperatures. Such a be- 
havior is inconsistent with the 1 /Ti ^ behavior 
expected from harmonic phonons at high tempera- 
tures due to the two-phonon Raman process [10]. 

In a previous work we succeeded in describing 
this anomalous temperature dependence consider- 
ing a strong anharmonicity of the vibration of the 
potassium ion [11]. We added a fourth order term 
to the harmonic potential and studied its infiuence 
within a self-consistent quasi-harmonic approxima- 
tion. Within this model the effective phonon fre- 
quency becomes a strongly increasing function of 
temperature, modifying the higher temperature be- 
havior of both NMR relaxation rate and resistivity 
[11], consistent with the experimental observations. 

In the present work we are studying a square- 
well potential as a model for the potassium dynam- 
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ics. This potential may seem less realistic and it 
also does not allow to "tune" the amount of anhar- 
monicity. However, it can be regarded as the most 
extreme case of "rattling" and, more importantly, 
the square-well potential allows us to calculate the 
phonon spectral function analytically, without hav- 
ing to resort to the quasi-harmonic approximation. 
We wish to demonstrate that an anomalous temper- 
ature dependence of the NMR relaxation rate like 
the one found in Ref. [11] results, reinforcing our 
previous conclusion. 

As the simplest possible model we consider the 
one-dimensional square-well potential with in- 
finitely high walls. As is well known, the energy 
levels in this case are given by 
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Here, L is the size of the well and m the mass of the 
atom. It is an easy exercise to analytically calculate 
all the matrix elements of the positional operator x, 
giving 
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if m -|- n odd and otherwise. The free, non- 
interacting phonon Green's function -Do('^) can be 
calculated quite generally from the expression [12] 
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Here, Z ~ q-Pe„ jg ^-j^g partition function and 
(3 = l/T the inverse temperature and we have 
normalized Do{u) to the lowest energy excitation 
(l|x|2). 

In the superconducting state the NMR relaxation 
rate at the potassium site exhibits a sudden de- 
crease. As has been pointed out in Ref. [9] this sug- 
gest that the rattling phonon must be strongly cou- 
pled to the conduction electrons. Following Ref. [11] 
we therefore introduce a finite phonon self energy 
n(a;) due to coupling to the conduction electrons. 
Then, the interacting phonon Green's function D{oj) 
can be calculated from Dyson's equation 
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The phonon spectral function A{lij) is obtained from 
1 
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Fig. 1. Phonon spectral function A(w) for the square well 
potential at three different temperatures T = 0.3£i (solid 
line), T = l.OBi (dashed hne), and T = 3.0Bi (dotted 
line) . With increasing temperature, spectral weight is shifted 
gradually away from the main peak at 3Ei to the higher 
energy peaks. 

Because the electronic energy scale is much larger 
than the phononic energy scale, we can make a low 
energy expansion of the phonon self energy and write 



n(a;) 
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where Hi and a are positive real constants describ- 
ing the interaction of the phonon with the conduc- 
tion electrons. Here, Hi is the renormalization of 
the phonon frequency and a determines the phonon 
damping rate. 

As can be seen from Eq. (3) the non-interacting 
phonon Green's function consists of a sum of poles 
at the transition energies £"„ — E„i ■ For a harmonic 
oscillator these poles all fall on top of each other, 
resulting in a single, temperature independent pole. 
However, for the present square-well potential the 
allowed transition energies are different and the cor- 
responding spectral function consists of a series of 
peaks at the positions Sn = 2mL'^ (^'^ "^^"^ 
ative weight of these peaks is temperature depen- 
dent. With increasing temperature spectral weight 
is gradually transfcrcd to higher energies, because 
higher energy transfers become possible. In Fig. 1 we 
show the energy dependence of the spectral function 
A{lo) calculated from Eqs. (l)-(5) at three different 
temperatures T = O.SE"! (solid line), T = l.QEx 
(dashed line), and T = 3.0£'i (dotted line). Here, 
El = is the ground state energy of the square 

well potential. For illustration, no phonon frequency 
renormalization Hi = and a small phonon damp- 
ing of a = 0.1 has been chosen. The influence of 
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these two parameters on the NMR relaxation rate 
will be discussed below. Fig. 1 illustrates how the 
spectral weight is shifted away from the lowest en- 
ergy peak at 3Ei to the higher energy peaks, when 

temperature is increasing. 

Within the quasi-harmonic approximation of 
Ref. [11] this gradual transfer of spectral weight to 

higher energies appears as a temperature increase 
of the effective phonon frequency. We wish to em- 
phasize that the peak positions do not change as a 
fimction of temperature. Only their intensities do. 
The effective phonon frequency of Ref. [11] there- 
fore corresponds to an effective average frequency 
and should not be associated with an energy level 
distance. For example, a neutron scattering experi- 
ment is expected to observe a temperature depen- 
dent intensity of a series of phonon peaks and will 
not be able to directly observe the effective average 
frequency. 

Once the phonon spectral function is known, the 
NMR relaxation rate 1 /Ti from quadrupolar inter- 
action of the potassium nucknis to the local elec- 
tric field gradient can be calculated. In Ref. [11] we 
have shown that the two-phonon Raman process is 
dominating over the direct process, because of the 
small number of phonon states available at the Lar- 
mor frequency [10]. Therefore, here we will focus on 
the two-phonon Raman process. In this case 1/Ti is 
given by 
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neglecting higher order vertex corrections due to 
electron-phonon coupling. Here, n{uj) is the Bose 
distribution function. In Fig. 2 we show a numeri- 
cal evaluation of Eq. (7) using the spectral function 
from Eq. (5). Again, IIi = and a = 0.1 has been 
chosen. It is immediately apparent that the NMR 
relaxation rate l/TiT possesses a peak as a function 
of temperature near T = 5Ei, in contrast to har- 
monic phonons. The behavior is qualitatively sim- 
ilar to the experimental data from Ref. [9]. At low 
temperatures \/T\T varies like due to the lin- 
ear behavior of the spectral function at low energies 
seen in Fig. 1. At high temperatures we can show 
that 1/T\T varies like 1/ s/T (see below). This is in 
contrast to harmonic phonons, where 1/TiT ~ T at 
high temperatures, but it is also difi^erent from the 
fourth order potential result where 1/TiT ~ const 
was found [11]. Apparently the high temperature be- 
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Fig. 2. Temperature dependence of the NMR relaxation 
rate 1/TiT from the two-phonon Raman process for the 
square-well potential. 

havior depends somewhat on the type of potential 
chosen. 

In order to derive the high temperature limiting 
behavior we assume that the damping parameter 
a is small enough such that the poles appearing 
in Eq. (3) remain well separated. Generally we can 
write Do{co) as a series of poles: 
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where a„ (T) are the temperature dependent weights 
of the poles. Because the matrix elements |(n|a;|m) |^ 
quickly decrease with increasing distance n — m, to 
a good approximation we can only keep the nearest 
neighbor terms m = n ± 1. Then we have 
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The partition function Z{T) = g-'^^i"' varies 
like ^/T at high temperatures because the number 
of terms contributing significantly to the sum in- 
creases as \/T- If the poles are well separated, the 
full phonon Green's function becomes 

OJ-Sn-an (T) n (w) 
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and the spectral function can be written in the form 
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with r„,e// = aan{T) (Sn - an(T)Ili). Using 
Eq. (11) and evaluating each pole separately, we 
find from Eq. (7) 
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where f2„ = 5n — a-a (T) Hi arc the renormalized 
phonon frequencies. Apparently, 1 /Ti scales like 1 /a 
for all temperatures and small values of a. This 
means that the parameter a has no large influence 
on the temperature dependence of the NMR relax- 
ation rate, but only on its absolute values. 

In the high temperature limit we find from Eq. (9) 

a„ (T) « ^^""l-^-^'^ ~ (13) 

Using n(f2„) ~ we finally arrive at 

— ^ ■ (14) 

We can estimate an absolute value of the expected 
peak position of the NMR relaxation rate by choos- 
ing the potassium mass for m in Eq. (1) and a size of 
the well of L « 1 A, as suggested by both bandstruc- 
ture calculations [1] and X-ray studies [2] . Then we 
find Ex « 6.1 K, which would suggest the NMR peak 
position to appear around 5i?i ~ 30 K. However, 
the peak in the experimental data appears near 14 K 
and also the peak is somewhat more sharp than the 
one in Fig. 2. A possible explanation for these dif- 
ferences could be the influence of a sizeable phonon 
renormalization due to electron-phonon interaction, 
as given by the parameter Hi in Eq. (6). To illus- 
trate the influence of the parameter Hi, in Fig. 3 
we show the temperature dependence of 1/TiT for 
four different values of IIi, keeping a = 0.1 fixed. 
With increasing IIi the peak becomes more sharp 
and the peak position moves to lower temperatures 
easily softening the peak position by more than a 
factor of 2. Thus, a better agreement with the NMR 
data on KOS2O6 could be obtained for a value of 
Hi w 0.5Ei. 

To summarize we have calculated the NMR re- 
laxation rate due to quadrupolar coupling to a lo- 
cal, strongly anharmonic phonon mode considering 
the two-phonon Raman process. We have calculated 
the free phonon Green's function analytically for a 
square-well potential. Including phonon renormal- 
ization due to electron-phonon interaction we were 
able to show that the temperature dependence of 
the NMR relaxation rate 1/TiT possesses a peak as 
a function of temperature, in stark contrast to har- 
monic phonons. The position and height of this peak 
depends on the strength of the phonon renormaliza- 
tion. We also derived the high and low temperature 




Fig. 3. Temperature dependence of the NMR relaxation 
rate 1/TiT from the two-phonon Raman process for different 
values of the real part of the phonon self-energy: IIi = Ei 
(solid Una), Hi = O.8B1 (dashed line), Hi = 0.5Bi (dotted 
line), and Hi = (dashed-dotted line). 

limiting behaviors. These results reinforce our pre- 
vious approximate treatment of a fourth order po- 
tential. 
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